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D . The method for determination the microscopic surface tension (t{R) of nanobubbles is developed 
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based on the new elaboration of the positronium bubble model. In contrast to existing structureless 
Ps bubble models, our version contains experimentally known molecular characteristics of liquids. 



■ The relationship, similar to Tolman's equation, between cr{R) and the radius R of the Ps bubble 
is derived on a microscopic basis. Numerical values for cr{R) are determined for a large number 
of liquids and liquified gases. The results are in agreement with the theoretical expectations and 
independent evaluations, available in literature. 
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I. INTRODUCTION 



More than one hundred years ago Willard Gibbs in the frameworks of thermodynamic approach 
estabhshed that if the surface of a hquid is curved, the surface tension coefficient a gets a function 
of curvature 1/r of the interphase boundary. Investigations of Gibbs were extended by Tolman 
Q and others It was found that behavior of oir) for droplets and bubbles is different. For 

droplets a{r) is approximately described by Tolman's equation: 

- T^r <^-^> 

where cxoo is the tension for a plane surface and the "Tolman distance" A is equal to the distance 
between the equimolecular dividing surface and the surface of tension. 

In case of bubbles o"(r) may pass through a maximum and than goes to zero, but in both 
cases at small r air) oc r ||^. Special additional assumptions are made concerning the behavior 
of the Tolman parameter A function of r fHl . 

However applicability of the thermodynamic approach to very small systems such as liquid 
droplets in the gas phase or bubbles in liquids remains questionable Basic equations hold for 
systems with large number of molecules. Besides it remains unclear how to relate locations of the 
dividing surface and the surface of tension to actual positions of molecules for nanosized voids and 
droplets. For this we need the microscopical theory of the surface tension. 

The problem is intricate because experimental verification of the theoretical assumptions made 
is practically impossible. Traditional techniques of surface tension measurements are not sensitive 
to the deviations of the microscopic surface tension at small r from the macroscopic one. Some- 
times the problem is masked by non-equilibrium properties of the surface formed, which depend 
on its "age" (so-called dynamical surface tension). 

Solvophobic property of the positronium (Ps) atom, i.e. its ability to form a nanobubble in 
liquid, makes it very attractive probe for investigation of the possible dependence of air) f^,^]. 
Recently we suggested new approach for determination of the microscopic surface tension in liquids 
using positron spectroscopy data Combining the data on the ortho-Ps lifetime and width of 
the "narrow" component of the AGAR (angular correlation of annihilation radiation) spectrum in 
the same liquid it is possible to extract parameters of the Ps bubble. Finally energy minimization 
condition allows to determine respective microscopic surface tension. 

In this paper we present detailed description of the Ps bubble model suggested in [§] and 
modification of the method of determination of the microscopic surface tension. 

Before starting the explanation of the subject we shell briefly remind the basis of the positron 



spectroscopy to introduce terms which will be used below. 



II. POSITRON ANNIHILATION IN MATTER 

The nonrelativistic approximation in QED gives the following expression for the spin-averaged 
probability per second for the annihilation of a positron-electron pair into two photons, for which 
the sum of the wave vectors k = ki + k2 is within d'^k 



dX = vrrgC ■ P27(k) 



d^k 
(2^ 



(2.1) 



Here tq = e'^/mc^ is the classical electron radius and c is the velocity of light. The photon-pair 
momentum density is 



P27(k) 



(2.2) 



where (/)+_(r_|_, r_) is the wave function of the annihilating e~^-e pair. To obtain total annihilation 
rate of the positron we should integrate Eq.(2.1) over the wave vector k, which gives 



/P27(k)^ = / \<P^4r,v)fd'r. 



(2vr) 



(2.3) 



In the case of intrinsic 27-annihilation of the Ps atom the wave function of the annihilating 
e~^-e~ pair may be written product 



^(r) ■ ^/'Ps(r^ 



(2.4) 



where ip{r) is the function of the center-of-mass coordinate r = (r_|_ + r_)/2 and ipps{r) depends 
on the relative coordinate r_|_ — r_. Fourier transform of Eq.(2.4) and integration over k in case 
when ippsir) = e'^^'^"'^ / ^Jstto^ is the ground state of the Ps (as = h/me^ is the Bohr radius) 
gives the following expressions for the photon-pair momentum density: 



P2,(k) = |#s(0)|^ 



P27(k) 



d^k 

(2^ 



l#s(0)p 



(2.5) 



From Eq.(2.5) and Eq.(2.1) we obtain the spin-averaged decay rate of the Ps atom A27 
It is equal to one fourth of the decay rate of the para-positronium Ap_Ps = rQc/2a^, Ap:';ps = 123 
ps.i 



^Account of the higher order corrections increases the p-Ps lifetime up to 125.2 ps. 



In case of e"*" annihilation in matter 0+_(r+,r_) should be replaced by a total positron and 
iV-electron wave function. In the frameworks of the independent particle model when all e+-e^ 
correlations are neglected can be approximately written as a product of the wave functions 
of the particles involved: 0+(r_,_) J2j 4'-\''^j)- Then P27(k) takes the form: 



(2.6) 



Here 0+ and are the unperturbed positron and electron wave functions, respectively, and the 
sum is taken over all occupied electron states. Integration of Eq.(2.6) over k gives 

Usually one assumes that 



P2,(k)^ = Ey ■ \e\r)\H'^r. (2.7) 



E ^ Z^^n, J |0+(r)|Vr = 1, (2.8) 

j 

where n is the number density and Z^s is the effective number of electrons per molecule capable to 
annihilate with the positron. Slow positron can not penetrate deep inside an atom, so core electrons 
do not contribute to Zcs- So Z^s is close to the number of the valence electrons. Substituting 
Eq.(2.8) into Eq.(2.7) and Eq.(2.1), we obtain Sommerfeld's result for annihilation rate of "free" 
positrons in matter 

K+ = nrlcZcsn. (2.9) 

III. PS BUBBLE MODEL 
A. Historical outlook 



In 1957 Ferrel [O suggested positronium in a liquid creates a nanocavity (bubble), repelling 



neighboring molecules outward. It happens because of a strong exchange repulsion between the 
electron, composing Ps atom, and electrons of host molecules. It was the onset of the Ps bubble 
model. Its two main aims are the calculation of the lifetime of the ortho-positronium (o-Ps) and 
calculation of the shape of the angular correlation of annihilation radiation (AGAR) spectrum, 
strictly speaking its "narrow component", which corresponds to intrinsic 27-annihilation of para- 
positronium (p-Ps). 

Residence of the p-Ps in a bubble practically does not change its lifetime. It is too short 
because of prompt intrinsic 27-annihilation. The situation is very different for the o-Ps state. 



Because of a restriction imposed by angular momentum conservation 27-decay mode is forbidden 
for o-Ps. Thus its most probable decay channel in vacuum is 37-annihilation. That is why o-Ps 
lifetime gets approximately 1000 times as large than of p-Ps one and constitutes 142 ns in vacuum. 
However in matter o-Ps may participate in another annihilation process. It is so-called pick-off 
process when e"*" annihilates into 27 with an electron of the opposite spin, belonging to surrounding 
molecules. Usually pick-off annihilation shortens o-Ps lifetime down to several nanoseconds. 

To account the pick-off annihilation process Tao modified Ferrel's model suggesting existence 
of an electronic layer inside the well, close to its boundary. Within the model, which uses an 
infinite spherically symmetrical potential well for simulation of the Ps bubble, introduction of 
such a layer is a unique way to account for the overlapping of the positron and external electrons, 
which is responsible for the pick-off process. Because of simplicity and physical transparency this 
model became very popular |ll2[ . 

By the end of the 50 's Stewart and Briscoe |T^ and Roellig |1^ introduced potential well of a 
finite height for more adequate simulation of the trapping potential of the Ps bubble. Their model 



came to present time practically without modifications. Its basis is the following [|T^. Ps atom 
is considered as a point quantum particle, which is self-trapped in a spherical free-volume cavity. 
Action of the surrounding molecules on the Ps is taken into account via an external potential, 
which is simulated by a spherical rectangular potential well with the depth U and radius Ru- 
A liquid is considered as a structureless continuum. Almost in all cases the energy of bubble 
formation is reduced to the surface energy, which is attributed to the interphase boundary. It is 
important that the position of this boundary is also associated with the location of potential well, 
i.e. with Ru. So the bubble formation energy is written as AnRfjaoo- In calculation of the pick-off 
annihilation rate it is assumed that external host electrons do not presented inside the potential 
well and o-Ps pick-off annihilation proceeds only due to overlapping of the o-Ps wave function 
with the electrons of a medium outside the boundary of the well. 



B. Smooth potentials and concept of average density 

Recently in |[T5|JT6|JT7|] the smooth potentials like f/ tanh^ (r/i?(/) and f/(l — e""^^^^) were used 
instead of the sharp finite well potential. With this the authors tried to take into account the 
smooth variation of the density of medium from the center of the bubble towards the bulk of the 
liquid. Application of these potentials is based on a possibility to obtain analytical expression for 
the Ps wave function and the energy of the Ps ground state. 

However we think that this concept of the "average" density profile in the problem of positro- 



nium formation is not well justified. One may admit that during the bubble formation stage the 
isolated molecules (so to say a "vapor phase") may exists inside a " pre-bubble" . But we can not 
take their presence into account in terms of smooth density distribution (and smooth potential). 
Ps motion is much faster then that of molecules and Ps wave function easily tunes up to cur- 
rent positions of molecules. Ps wave function goes to zero inside the molecules because of strong 
exchange and coulombic repulsion. Such a behavior of the wave function increases Ps kinetic 
energy (zero-point energy) and finally promotes pushing all the molecules out of the cavity to 
the boundary of the Ps bubble. In contrast to empty bubbles (no Ps inside), the equilibrium Ps 
bubble does not contain the "vapor phase". Presence of the light quantum particle in the bubble 
leads to the much more abrupt density profile on the boundary of the Ps bubble. Only in this 
case the usage of the bell-like Ps wave function (similar to that in Eq.(3.1)) is meaningful. We 
conclude that structure of the interphase region of the Ps bubble is different from that of usual 
vapor-liquid boundary. 

Obviously smooth Ps wave function is not good approximation outside the bubble. Self- 
consistent results should correspond to a rather small penetration of the Ps wave function to the 
bulk of a liquid. As we shell see below the results obtained on the basis of the present model are 
well- matched with this requirement. 

Of course the Woods-Saxon potential is suitable for smoothing the sharp edge of the potential 



of the rectangular wall \TE\. In the most realistic case a <^ Ru (a is the third parameter (beyond 
Ru and U) entering the Woods-Saxon function) the potential approaches to the square well shape 
and the results of the fitting of experimental data using these two potentials should be similar. 



C. General comments about modifications of the Ps bubble model 

Basing on the above comments, we adopt here the potential of the finite spherically- 
symmetrical rectangular well for simulation of the Ps bubble, but introduce additional specifi- 
cations to make the bubble model more realistic. 

As we have seen in the standard model Ru is overloaded by different physical mean- 

ings. It determines the position of the potential well, which confines Ps in the bubble. At the 
same time Ru determines the bubble formation energy, A-KR^a^o- Ru is also related to the pick-off 
annihilation rate of Ps. It is clear that description of these effects having such a different physical 
nature by means of one parameter is very crude. So we attempt to split theses effects, introducing 
additional parameters. Below we reserve for Ru the meaning of the position of the potential wall 
only (Fig.l). 



How do calculate the energy of the bubble formation? It is an important question for the Ps 
bubble model. In section 5 we shell see that the naive estimation AtiR^jCToo for the surface energy 
contribution is not correct. Position of the potential wall responsible for "reflection" of the Ps 
into the bubble and position of the surface of tension, related to intermolecular interaction, are 
different. They are placed on different distances from the center of the bubble. Approximating 
molecules by spheres interacting with each other by means of "central" forces, it seems reasonable 
to identify the surface of tension with the sphere Srj^r^^ passing through the centers of molecules 
residing on the first molecular layer of the interphase boundary (Fig.l). Here R is the radius of 
a free- volume and -Rws is the Wigner-Seitz radius (|7ri?^g = l/n).l In typical Ps bubbles the 
difference between R and R + R^is is important for estimation of the surface energy. This problem 
will be discussed in Section 4 in more details. 

To calculate the o-Ps lifetime we need to know distribution of the electronic density close to the 
boundary of the Ps bubble. Depending on U and Ru, may penetrate in some extent through 
the nearest molecules (bubble boundary) into the bulk of the liquid. At the same time electrons 
belonging to the nearest molecules may reside inside the potential well. Thus, electronic density 
profile and that of the potential, localizing Ps, should not coincide. Our approach explicitly takes 
into to account penetration of outer electrons into the bubble through the parameter 5 (Fig. 1). 
It leads to additional pick-off annihilation within the layer of the thickness 5 close to Sij^-sphere, 
but inside it.i It is worth noting that this circumstance allows to reproduce standard infinite 
potential well bubble model as a limiting case oo of our approach. 

D. ortho-Ps lifetime 

By the end of bubble formation (when all the molecules of the "vapor phase" are pushed out 
to the bubble boundary) the Ps center-of-mass wave function takes the form 



^If molecules are not spherical, but rather elongated, it could be reasonable to approximate them as a 
equence of spherical fragments. In this case |'?i"-Rws takes sense of the volume of the fragment. 

^Of course there are some other reasons, which could lead to deviation between density and potential 



profiles. For example, approaching to the boundary of the bubble one may expect small deepening of 



the potential due to polarization interaction as was discussed by Chuang and Tao |19| in application to 
silicagel powders. We neglect such effects here. 



sin kur 



il){r) 



H 



r <Rj 



27r(l + kR, 



u) 



r 

(3.1) 



where 



kfj = 2mp,E/h^ = 2mps(f/ - E) /h^ = 2mp,U/h^ - kfj. (3.2) 

Here E is the kinetic energy of Ps and mps = 2me is its mass. Requirement of smoothness of iplr) 
at r = Ru leads to 

K = —kjj cotkuRu, n/2 < kuRu < TT. (3.3) 

Energy spectrum of the Ps in the well can be obtained form this equation. When kuRu ^ vr/2 
and >iRu the Ps ground state escapes from the potential well, while the limit kuRu vr and 
>cRu —>■ oo corresponds to the infinite potential well. The relationship 

U = , ^ (3.4) 

Zmps sm kjjRu 

follows from Eq.(3.2) and Eq.(3.3) and will be used below. 

The rate of the pick-off annihilation can be roughly estimated approximating the factor 
Y,j |0^^(r)p in Eq.(2.7) as ^efr'^ ■ ^{^ > R)- Here -i^-function equals to unity, if r > i?, other- 
wise it is 0. Than the r.h.s. of Eq.(2.7) reduces to 

^ J \Mr)\' ■ |0^'^(r)r ■ d^r ^ Z^suPr, Pr = mr)f d^r. (3.5) 

Here we approximated 0+(r) by the Ps center-of-mass wave function iplr). Pr is a probability to 
find Ps (and therefore e+) outside the free-volume sphere Sr. Thus we obtain the relationship for 
the pick-off annihilation rate of the o-Ps atom §: 

Ap_ofT = TcrlcZ^snPR. 

In small bubbles Api^^ practically coincides with the o-Ps lifetime To_ps, but in large bubbles (for 
example in liquid He) we should take into account intrinsic 37 decay of the o-Ps, which proceeds 
with the rate Aa^ = 1/142 ns^"*^: 



^Strictly speaking the annihilation rate of e"'"-e~-pair having zero spin is 4 times as large, but the number 
of host electrons which may form such a zero-spin pair is 4 times less. Thus, these effects cancel each 
other. 



It is convenient to separate Pr it into two parts: 



(3.6) 



Pr = Ps + P 



'ru, Ps= f mr)fd^r, Pr^ = T \i;{rf d^r, Ru = R + 6. (3.7) 

'J R 'J Rjj 



rR+S 

l^l^irM'd-'r, Pn„ = , 

IRu 

Straightforward integrations and account of the normahzation condition for il){r) give 

p _ kuS - sin ku6 ■ cos{2kuRu - kyS) ^ _ sin^ kyRu 

kuRu — tan kuRu ' ^ 1 — kuRu cot kuRu 

In the limit of the infinite potential well {kuRjj vr, xRjj oo, Pr^^ 0) Eq.(3.8) for Ps is 
reduced to the well-known Tao formula for the o-Ps lifetime |jl2| 



= = A _ !!!i(MZM, (3.9) 

To-Ps Ru 27r 

where r^.pg is usually identified with the positronium spin-averaged lifetime 0.5 ns. The infinite 
potential well approach is very popular because of its simplicity. Knowing experimental value 
of the o-Ps lifetime and calculating Ru from the AGAR data (Eq.(3.21)), Eq.(3.9) may give an 
information about 6. 

In the finite potential well model we suggest to define parameter S in the following way. It was 
noted by Kobayashi |2^, that if in Eq.(3.9) the free- volume radius R = Ru — S tends to zero and 



6 = 1.66 A, the energy of the Ps ground state gets equal to the Ps binding energy in vacuum (6.8 
eV) . It indicates that Ps may not exist without the free volume. Generalization of this hypothesis 
for the case of the finite well is the following: when i? — or Ru — > S, the Ps bound state escapes 
from the potential well. It leads to the following relationship between U and 6: 

u.E{,.u.-..s-l)-^^.'-lny{^)\ (3.10) 

where Ry= = 13.6 eV. Of course in an unperturbed liquid (without the bubble) e"*" and e 
remain bound due to the long-rage Goulombic interaction, but their binding energy will be small 
and the average e"'"-e~ distance becomes larger than intermolecular distance. It is the quasi-free 
positronium in matter. Sometimes it is called as the swollen Ps. Electron density of the "own" 
electron on the positron in such a state is small in comparison with the other electrons. So, 
positron annihilation will look like the free e"*" annihilation. 
Gombination of Eq.(3.10) and Eq.(3.4) gives 



kuS = — sin kuRu- (3-11) 

This relationship essentially simplifies our approach, because now becomes a function of the 
parameter kjjRu only. Thus, knowing o-Ps lifetime we can directly obtain kjjRjj. However, 
an additional information is needed to obtain Ru and ku separately. For this purpose ACAR- 
spectroscopy data will be used. 



E. Narrow component of AGAR spectra 

The distribution of annihilating photons over k^ is measured by means of the long-slit angu- 
lar correlation annihilation apparatus. The most reliable information about the Ps state in 
the bubble is obtained from the shape of the narrow component of the AGAR spectra, which 
corresponds to intrinsic annihilation of the p-Ps: 

/+00 r+oo 
dk^ / dkyP2^i}i). (3.12) 

Calculation of the Fourier transform of tp{r) Eq.(3.1) gives the photon-pair momentum density: 



P2-y{k) OC 



1.2 _ U2 



sin kuRjj cos kRi 



k 



COS kuRu sin kRu + 



If X 

— sin kuRu cos kRu + — sin kuRu 

k^ + \ k 



(3.13) 



For the infinite potential well it reduces to 



p2^{k) OC 



kfj siv? kRu 



(3.14) 



F {kl-k^f 

It is convenient to carry out an integration over kx- and fc^-components of the photon wave vector 
using the following transformation, k"^ = k1 + k^: 



+00 



dkrf 



+ 00 



dkiL 



/■oo roc r 

27r / kx.dk^ . . . = vr / dik^ - kl) . . . = t\ 

Jo Jo Jk 



Integrating expression for p2-f{k) in such a manner, we obtain ^1 



oo 

TT / kdk . 



(3.15) 



roo 

N{kz) OC / kp2^(k)dk oc 
which in the limit of the infinite well gives 



P cos P + hRu sin [3 

{kiRi - p^){Kmi + p^ 



(3.16) 



d(3 sin^ [3 



(3.17) 



The value of the full width at half maximum, Ofwhm? of the narrow component of AGAR 
spectrum (see Eq.(3.16)) can be obtained from the following integral equation 2iV(A;FvvHM/2) = 
N{0) or 



2 / ^{(3, kuRu)d(3 = / ^{(3, kuRu)df3, 
Je, Jo 



^{13, x) = - 



where 



^FWHM^ 



(3 cos P — xsinP cot x 
{x^ - (3^){(3^ + cot^x) 



nT-eC Ofwhm 

-Ku- 



(3.18) 



(3.19) 



2 h 2 

Here mgC is a momentum of one of the annihilating 7-quanta. 

It is important that knowing kjjRu (from the o-Ps lifetime data) and 0fwhm from AGAR 
measurements we may obtain Ru and all other parameters of the Ps trap {U and 5). If [/ — 00, 
in Eq.(3.18) we should set x — * vr, which leads to 



°° d(3 ( sin 13 



e. (3 \7T^-(3^ 



dp ( sin P 



P \7T^-P^ 



(3.20) 



with the same meaning of 6*. Numerical solution of Eq.(3.20) gives 6* = 2.1480. Substituting 
this value to Eq.(3.19) we obtain simple relation between 6fwhm and radius of the bubble 



Roo [A] 



16.65 



(3.21) 



©FWHM [mrad]' 

Usually for extraction of the narrow component from the total AGAR spectrum it is decom- 
posed into a set of gaussians. To visualize uncertainty which comes from neglecting deviation 
in shape between Eq.(3.16) and the gaussian with the same Bfwhm, in Fig. 2 we plotted several 
spectra for different Ps traps. For rather "deep" well {kuRu > 2.5) the difference is not large, but 
for "shallow" traps {kuRu < 2.5) this deviation should be taken in to account in decomposition 
of the AGAR spectrum. 



IV. ELEMENTARY MODEL OF A CAVITY FORMATION 



Let molecules in the liquid interact, for example, through to the Lennard- Jones potential: 



When r is equal to average molecular separation f the pair- wise molecular potential energy reaches 
its minimum: 93Lj(r = f) = — e. In equilibrium f is related with the number density n as 
f fa We shall consider molecules as semi-hard spheres of radius f/2 with only nearest 

neighbor interactions {^Plj approaches zero very rapidly with increasing r). 

Since the energy required to form a surface arises from the decrease in number of bonds of 
molecules at the boundary, the accompanying decrement in coordination number, z/, needs to be 
estimated. The bulk liquid we imagine as a closely packed structure. Everywhere in it one may 
find four nearest molecules, which form a regular tetrahedron of edge length f, whose centers lie on 
a sphere of radius ri = t^?,/'& (Fig. 3). Therefore, the area per molecule on Sr^ is \ -A.T:r\ — irrf 
and the coordination number will be 

In what follows we neglect the difference between ri and the Wigner-Seitz radius because it may 

1 /3 

be seen that ri = (^^^ ■ -y/|i?ws = 0.986i?ws- Hence for all practical purposes we may take 
that the volume of the sphere of the closest neighbors to have the same volume as the average 
volume per molecule in the liquid or, in other words, ri ~ Rws- 

Next permit the positronium to create a spherical free volume AttR^/S (or bubble) of radius 
R inside the tetrahedron and thereby pushing the molecules outward. The other molecules also 
rearrange themselves to settle on the first molecular layer (FML) of radius ri + R, Fig. 4. The 
number A^fml of molecules lying on this layer can be estimated from its area, 47r(ri-|-i?)^, dividing 
by the area per molecule, irrf, that is: 

Sri+R . fri + R^ ^ 



iVFML = ^ = 4 . (4.3) 

Si \ ri / 

Because of formation of the bubble with the free volume AttR^/S the coordination number which 
was uq suffers a decrement in proportion to the free area per molecule residing on the FML, 
Sr/Nfml: divided by the area per molecule in the bulk, 47r(f/2)^, namely 



A. = ^^^^-o=4(^) , Sn = AnR\ (4.4) 
ttH Vri -I- R 



The fraction f^l^2)^ represents a factor decreasing the coordination number. Accordingly, the 
surface energy should be proportional to the number of broken bonds viz. 



Ea ~ NfmlAu ~ R^/rl 



(4.5) 



Also in view of the fact that we are deahng with the central forces and these act at the center of 
molecules we expect the surface of tension is located aX R^j = R + ri ^ R + Rws- Thus one would 
expect the surface energy to be proportional to the area of the surface of tension, Sr^ — 47ri?^, 
with the curvature dependent coefficient of proportionality a{Ra-)' 

E,r^Rla{R„). (4.6) 

In the limit — > oo the surface tension energy E^j should reproduce standard relationship for 
the plane surface viz. E^ — AnR'^a^. So we reconstruct coefficient of proportionality in Eq.(4.6). 
Finally, comparing Eq.(4.5) and Eq.(4.6) we obtain 



E, = A7rRla{R,), ^ (1-)' ^ J- R, = R + n ^ R + R^s- (4.7) 

(Too V-rio-/ (i + itws/-n'j 

Furthermore, regarding subsequent minimization of the total energy (balance condition), we 
prefer to rewrite the surface energy of the bubble in an integral representation for convenience 
and thus introduce the surface tension function a{r) through: 



r 

which also happens to be the Laplace form. Of course if a is put equal to and integration is 
performed over the volume AnR^/S, one obtains the result E^- — At^R^Goo- However, here we must 
proceed with the space integration over the range from i?ws to R^ and obtain d{r) by solving the 
integral equation: 



d-^r = 47ri?Voo. (4.8) 

'iiwS f 

Eq.(4.8) being differentiated with respect to i?, gives us a Tolman-like expression: 

This relationship sheds light on a physical meaning of the Tolman length A in Eq.(l.l) through 
the relation 2 A ^ Rws- Thus from the molecular point of view A accounts for the molecules on 
the curved first molecular layer as having more neighbors (less broken bonds) that the molecules 
on the plane interphase surface. It is worth noting that the surface tension coefficient of the 
curved boundary depends not only on the position of the surface of tension, but also on the type 
of representation of the surface energy. 



V. ENERGY MINIMIZATION 



Energy minimization condition can be naturally inscribed into our consideration in the fol- 
lowing manner. As we have demonstrated above all parameters of the Ps trap can be obtained 
from the o-Ps lifetime and the width of the narrow component of the AGAR spectrum. However 
a, entering Eq.(4.8), may be considered as an unknown function of R, neglecting the theoretical 
prediction, Eq.(4.9). Thus we suggest to extract so-to-say "experimental" dependence of 5"(i?) 
from the principal of the minimum of the total energy Etot of the Ps bubble!: 



tot 



E, + ^R'p^ + E= I 
3 Jh 



3 2mps 



(5.i: 



Here we added the term ^R^Pq, which is the work against external pressure Pq. This term is 
important for rather large Ps bubbles in liquified gases. 

To proceed with the minimization of Etot let us first figure out a useful relationship for dE/dR: 



dE 



2kuE 



(5.2) 



dR tan kuRu ~ kjjRu 
which can be obtained from Eq.(2.2) and Eq.(2.3) by differentiation over R and further exclusion 
of X. Then the balance condition dEtot/dR = may be written as 
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This equation gives 
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For the infinite potential well Eq.(5.4) simplifies to: 
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(5.5) 



The usage of the integral representation for allows to avoid an appearance of the derivative of a 
over R minimizing Etot- 



VI. RESULTS AND DISCUSSION 



A. The bubble is a deep potential well for Ps 

Using Eqs. (5.4-5.5) and knowing Tq-ps and Ofwhm from the e~^-e~ annihilation experiments it 
is possible to obtain microscopic values of a for the interphase boundary having curvature radius 
about several angstroms.i Deriving these relationships we did not use any particular expression 
for the surface tension coefficient vs curvature radius of the boundary. We assumed only that 
cavity formation energy is written in the Laplace form, Eq.(4.8). 

For all investigated molecular liquids ratios o joaa are less than unity. Variations between the 
values corresponding to the finite well, Eq.(5.4), and the infinite well, Eq.(5.5), is small (Tables 1, 
2; Fig. 5). The reason is that the obtained depth V of the potential well is rather large U ^ E. In 
Section 3 we mentioned that this inequality ensures self-consistent usage of the bell-like Ps wave 
function in the frameworks of the present theory. 

In [1^ in some liquids (glycerin, ethylene glycol, methanol- water mixtures) was not fulfilled. 



It implies large penetration of the Ps to the bulk of a liquid, which means that the results might 
not be reliable.! 



B. Profiles of the potential well and that of electronic density are different 

Previous formulations of the Ps bubble model, utilized finite well potential, were not able to 
reproduce in a limiting case U ^ oo the Tao formula (3.9) for Tq-ps, obtained within the infinite 
potential well model.i We avoid this drawback in the present formulation through the parameter 
6, which accounts some penetration of the outer electrons within the well. Thus 6 discriminates 
profiles of the potential and that of electronic density. 



^Without the last term with pQ Eq.(5.5) was obtained in 

''We think that in |jl5,16,17,lS] expression for the pick-off annihilation rate the square of the Ps psi- 
function | "0(^)1^ has to be multiplied on the respective value of the electronic density and than this 
product should be integrated over space variables. 

^In the frameworks of the conventional finite well model in the limit C/ — > oo the Ps wave function 
is confined within the bubble only and does not overlap with outer electrons. Therefore the pick-off 
annihilation rate equals to zero. 



C. Separation of the position of the potential wall and the surface of tension 



Another new and important element of the present formulation is the separation of the position 
of the potential well (Ru), which reflects Ps into the bubble, and that of surface of tension 
{R + Rws): related to interaction!! between the molecules residing on the first molecular layer. 

D. Correlation between the Tolman length and the Wigner-Seitz radius 

Comparison between "experimental" values of ojo^a and respective theoretical prediction, 
Eq.(5.9), for two values of the parameter 2 A (-Rws and 3-Rws) are shown in Fig. 5. One may 
conclude that there is reasonable agreement between the theory and experimental data in spite 
of many simplifications done, which could be inadequate especially in polar liquids or in liquids 
with large non-spherical molecules. 

Assuming validity of the Tolman equation and knowing d ja^ values, one may calculate the 
ratios 2A/R-ws. They are within the interval from 1 to 3, which is in a reasonable agreement with 
Eq.(4.9).0 Our values of 2 A for water are 5.4 and 3.3 A (see Tables 1 and 2), which correlates 
well with available literature data0: 6.0 p2[, 1.8 and 2.0 ^ A. The same takes place in liquid 



argon: we obtained 2A = 4.5-3.8 A, while in 2A = 7.3 A. 

In spite of the large uncertainty of the A/R-yjs values it is worse noting their correlation within 
the classes of different chemical compounds. In isooctane, neopentane and tetramethylsilane, i.e. 
in liquids with round molecules values of 2A/i?ws are close to unityEl in agreement with Eq.(4.9). 
In liquid hydrocarbons made up from normal, cyclic and aromatic molecules, in higher alcohols, 
diethylether, acetone the ratio 2A/i?ws increases up to 1.8 in average. Probably it is related to the 
orientation of molecules when their maximal linear dimension primarily directed perpendicularly 
to the surface of the Ps bubble. Higher values of 2A/i?ws (up to 2.7) occur in low alcohols, water 
and acetonitrile. These liquids consist of small polar molecules strongly interacting with each 
other. It is interesting that 2A/_Rws in CS2 is also high in spite of CS2 molecule has no dipole 
moment. However, as follows from radio-spectroscopy studies [EHl, complicate polar molecular 



^rupture of intermolecular bonds. 

^'^Experimental uncertainty of 2A/i?ws is about 100%. 
^^Data for droplets. 

^^Small value of 2A/i?ws in 1,4-dioxane, which does not belong to this class of compounds, probably 
related to the special alignment of dioxane molecules on the surface of the bubble. 



associates are presented in liquid CS2. As a result, an effective dipole moment per molecule turns 
out to be approximately 0.1 D and intermolecular binding energy between CS2 gets about 0.05 eV. 
Thus one may expect that there is a correlation between 2A/i?ws and efficiency of intermolecular 
interaction. The latter can be characterized, for example, by critical pressure p^r — i^-, where 
a and h are known parameters of the Van-der-Waals equation (Fig.6). This correlation is clearly 
seen in a homological series of alcohols. It is difficult to expect better correlation between the data 
obtained from rather schematic Ps bubble model and parameters of the Van-der-Waals equation, 
applied to the liquid phase. 

E. Ps is an electrically neutral probe of the interphase boundary 

A perturbation of the surface caused by the presence of the Ps atom in the bubble on does 
not extend deep inside the liquid because of electrical neutrality of the positronium. One may 
expect that the structure of the interphase boundary in the Ps bubble will be more close to the 
free surface than in the case of the bubbles formed by excess electrons. It is also worse mentioning 
that the surface of the Ps bubble is rather "fresh". Its age is no more than some nanoseconds. 
Contrary, surfaces studied by means of conventional methods have the ages many orders higher. 
So in liquids with rather long relaxation times properties of the boundary of the Ps bubble and 
that of the equilibrium surface may be different. 

VII. CONCLUSION 

Major part of this work is the development of the Ps bubble model, which forms the basis of the 
method for determination if the local surface tension. The modifications done are the following. 

Firstly, it is taken into account, that position Ru of the potential well docs not coincide with the 
position of the surface of tension. Secondly, we admit a possibility of the electrons of the nearest 
molecules to penetrate inside the Ps bubble. Just this feature allows us a to reproduce the Tao 
formula as a limiting case C/ — > 00 of our finite potential well model. It is important that we did 
not introduce undefined parameters to the model. It is due to the additional constrain, Eq.(3.10), 
which has the following physical meaning: in molecular liquids "preexisting" free volume can not 
localize Ps atom; in an "unperturbed" liquid Ps exists in the quasi-free (swollen) state which 
manifests experimentally like free e"*" annihilation. 

Elementary "geometric" consideration of the cavity formation, done on a molecular level, 
made possible to reproduce the Tolman equation and clear up the physical sense of the curvature 
dependence of surface tension. We have found that Tolman's length A takes into account an 



increase of the number of the nearest neighbors of a molecule, residing on the first molecular 
layer of the curved boundary of the bubble in comparison with coordination number of a molecule 
on flat interphase boundary. It is shown that particular value of the surface tension coefficient 
depends on the position of the surface of tension, and on the choice of explicit expression for the 
energy of the bubble formation. 

It is shown that the data on o-Ps lifetimes and widths of narrow component of AGAR spectra 
allow to obtain "experimental" values of the surface tension coefficient without any hypotheses 
about its concrete functional dependence vs R. The results are in a satisfactory agreement with the 
Tolman relationship Eq.(4.9) and other independent evaluations of the surface tension coefficients. 
We have found a correlation between the ratio 2A/i?ws a-nd critical pressure. It indicates on its 
usefulness in consideration of intermolecular interactions and structure of surface layers. 
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Figure Captions 



Figure 1. 

The Ps bubble. The ccntcr-of-mass wave function ip{r) of Ps confined by a spherical potential 
well of the depth U and radius Ru, Eq.(3.1). AttR^/S is the free volume. R + i?ws is the distance 
from the center of the bubble to the centers of molecules residing on the first molecular layer. The 
parameter 5 characterizes penetration of host electrons inside the potential well of the bubble. 

Figure 2. 

Different normalized (unit area below lines) narrow components of AGAR spectra, having the 
same ©fwhm- Circles represent the gaussian line ^^^^^^^ exp(— 2^), © = hkz/meC with the second 
moment M2 = ©|^jjM/(81n2). Solid line represents the narrow component for in the case of 
an infinite potential well, Eq.(3.17). Dashed hues are plotted according to Eq.(3.16) for different 
values of kuRu- 

Figure 3. 

Tetrahedron formed by joining the centers of four nearest neighbors which in turn lie on a sphere 
of radius ri (the fourth molecule, nearest to the reader, is not shown). Molecules are simulated 
as spheres of radius f/2 (f being the average intermolecular distance). 

Figure 4. 

Creation of the free volume spherical void in a liquid. 
Figure 5. 

Dependence of the relative surface tension vs i?/it!ws- Values of a/a^ for different liquids are 
represented with a help of the respective numbers listed in Table 2. Upper line represents Tolman's 
equation (4.9) and the curve below is the same relationship but with the factor of three in the 
Tolman parameter viz. a/aoo = (1 + 3i?ws/-R)~^- 

Figure 6. 

Correlation between 2A/i?ws a-nd critical pressure Pcr in different liquids at room temperature 
(Table 2). Linear proportionality between A/i?ws and Pcr in homological series of alcohols is 

clearly seen (corresponding numbers are encircled). 
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TABLE I. Some parameters of liquids and Ps bubbles, obtained within the infinite potential well 
model. 



liquid 


A 


em 2 


'0— Ps 
ns 


fc'FWHM 

mrad 


7? 

-'•'OO 

A 


A 

u 

A 


A 


a 


2A 

Rws 


n-C5Hi2 


n-pentane 


3.583 


15.32 


4.25 


2.25 


7.4 


2.0 


5.4 


0.43 


2.0 




n-hexane 


3.737 


17.74 


3.92 


2.26 


7.4 


2.1 


5.3 


0.38 


2.3 


n-C7Hig 


n- heptane 


3.881 


19.65 


3.85 


2.33 


7.1 


2.0 


5.1 


0.37 


2.2 


n-CioH22 


n-decane 


4.267 


23.3 


3.49 


2.47 


6.7 


2.0 


4.7 


0.38 


1.8 


n-Ci2H26 


n-dodecane 


4.492 


24.84 


3.43 


2.43 


6.8 


2.0 


4.8 


0.33 


2.2 


n-CuHao 


n-tetradecane 


4.696 


26.0 


3.35 


2.52 


6.6 


2.0 


4.6 


0.35 


1.9 


i-CgHis 


isooctane 


4.038 


18.33 


4.05 


2.42 


6.9 


1.9 


5.0 


0.45 


1.5 


C(CH3)4 


neopentane 


3.656 


11.52 


5.15 


2.21 


7.5 


1.9 


5.6 


0.53 


1.4 


C6H12 


cyclohexane 


3.509 


24.65 


3.24 


2.45 


6.8 


2.1 


4.7 


0.38 


2.2 


C7H14 


methylcyclohexane 


3.692 


23.28 


3.50 


2.49 


6.7 


2.0 


4.7 


0.41 


1.8 


CeHg 


benzene 


3.285 


28.22 


3.15 


2.55 


6.5 


2.0 


4.5 


0.39 


2.1 


C6H5CH3 


toluene 


3.486 


27.92 


3.24 


2.57 


6.5 


2.0 


4.5 


0.40 


2.0 


C2H5C6H5 


ethylbenzene 


3.654 


28.74 


3.02 


2.44 


6.8 


2.1 


4.7 


0.32 


2.8 


(CH3)2C6H4 


o-xylene 


3.769 


29.76 


3.08 


2.54 


6.5 


2.0 


4.5 


0.35 


2.3 


(CH3)2C6H4 


m-xylene 


3.658 


28.47 


3.20 


2.49 


6.7 


2.0 


4.7 


0.34 


2.4 


(CH3)2C6H4 


p-xylene 


3.663 


28.01 


3.21 


2.49 


6.7 


2.0 


4.7 


0.35 


2.4 


(CH3)3C6H3 


mesitylene 


3.748 


27.55 


3.21 


2.49 


6.7 


2.0 


4.7 


0.35 


2.3 


(CH3)4C6H2 


1,2,3,4-thetra- 
methylbenzene 


3.888 


29. 


3.02 


2.52 


6.6 


2.0 


4.6 


0.34 


2.2 


CeFe 


hexafluorobenzene 


3.572 


22.63 


3.78 


2.39 


7.0 


2.0 


5.0 


0.37 


2.4 


Si(CH3)4 


tetramethylsilane 


3.780 


13.20 


4.75 


2.25 


7.4 


1.9 


5.5 


0.49 


1.5 


(C2H5)20 


diethylether 


3.453 


16.65 


3.82 


2.29 


7.6 


2.1 


5.2 


0.44 


1.9 


1,4-C4H802 


dioxane 


3.234 


32.61 


3.02 


2.86 


5.8 


1.8 


4.0 


0.52 


1.1 


CH30H 


methanol 


2.528 


22.12 


3.58 


2.29 


7.2 


2.1 


5.2 


0.37 


3.5 


C2H50H 


ethanol 


2.855 


21.97 


3.50 


2.35 


7.1 


2.1 


5.0 


0.39 


2.7 


C3H70H 


propanol 


3.101 


23.32 


3.38 


2.40 


6.9 


2.0 


4.9 


0.39 


2.5 


C4H90H 


butanol 


3.311 


24.93 


3.36 


2.46 


6.7 


2.0 


4.7 


0.39 


2.3 


CgHi/OH 


octanol 


3.938 


27.10 


3.13 


2.57 


6.5 


2.0 


4.5 


0.39 


1.8 


H2O 


water 


1.928 


72.14 


1.85 


3.05 


5.4 


2.0 


3.4 


0.39 


2.8 


D2O 


heavy water 


1.930 


70.89 


1.95 


2.87 


5.8 


2.1 


3.7 


0.31 


4.1 


(CH3)2CO 


acetone 


3.078 


24.02 


3.29 


2.45 


6.8 


2.0 


4.8 


0.41 


2.2 


CH3CN 


acetonitrile 


2.756 


28.66 


3.30 


2.44 


6.8 


2.0 


4.8 


0.35 


3.2 


CS2 


carbon disulfide 


2.880 


31.58 


2.20 


2.35 


7.1 


2.5 


4.6 


0.29 


3.9 


He, 4.2 K 


helium 


2.350 


0.096 


99.1 


0.86(6) 


19.3 


1.9 


17.4 




0.2tl 


H2, 20.3 K 


hydrogen 


2.243 


1.92 


28.6 


1.3(2) 


12.8 


2.1 


10.7 


n 44+0-34 


6.1^1" 


N2, 77.3 K 


nitrogen 


2.395 


8.85 


11.0 


1.8(2) 


9.2 


1.8 


7.4 






Ar, 86.4 K 


argon 


2.245 


12.42 


6.50 


2.20(15) 


7.6 


1.8 


5.8 


o.56l°:}i 


2.0tl 



i?ws is the radius of the Wigner-Seitz cell (calculated from the density and molecular mass). 



is the surface tension coefficient of a liquid with plane interphase boundary at room tempera- 
ture (except the cases of liquified gases). 



To-Ps is the o-Ps lifetime [12,21 . 



©FWHM is the full width at half maximum of the narrow component of the AGAR spectra [0,0]. 
i?oo = R + S is the radius of the Ps bubble, calculated from Eq.(3.21) using 6fwhm values. 



S is the penetration depth of the outer electrons into the Ps bubble. It is obtained from Eq.(3.9) 
using Roo and Tq-ps values. T°_pg was adopted to be 0.5 ns in all cases except He (1.9 ns), H2 
(0.92 ns) and N2 (0.56 ns). 

Relative uncertainties of To_ps and 0fwhm are no more than 5 and 10% respectively (for liquified 
gases indicated in parenthesis). Uncertainty of a /aoo is approximately four times larger than un- 
certainties of 0FWHM (about 40%). 

2A/i?ws are calculated from Eq.(l.l) using respective values of a/cxoo and assuming r = i? (in 
accord with Eq.(4.9)). 



TABLE II. Parameters of the Ps trap and microscopic surface tension obtained within the finite 
potential well model. 





liquid 


Ru 
A 


5 
A 


R 
A 


U 
eV 


E 
eV 


kuRu 


a 


dyn 
cm 2 


2A 
Rws 


Per 

MPa 


1 


n-pentane 


6.4 


1.2 


5.2 


3.48 


0.37 


2.81 


0.49 


7.45 


1.5 


3.36 


2 


n-hexane 


6.3 


1.2 


5.1 


3.32 


0.38 


2.80 


0.43 


7.56 


1.8 


3.01 


3 


n-heptane 


6.1 


1.2 


5.0 


3.49 


0.40 


2.80 


0.42 


8.31 


1.8 


2.76 


4 


n-decane 


5.7 


1.1 


4.6 


3.66 


0.45 


2.78 


0.43 


9.99 


1.4 


2.10 


5 


n-dodecane 


5.8 


1.2 


4.6 


3.52 


0.44 


2.78 


0.37 


9.24 


1.8 


1.83 


6 


n-tetradecane 


5.6 


1.1 


4.5 


3.72 


0.47 


2.78 


0.40 


10.3 


1.4 


1.61 


7 


isooctane 


5.9 


1.1 


4.8 


3.91 


0.43 


2.80 


0.51 


9.27 


1.1 


2.57 


8 


neopentane 


6.6 


1.1 


5.5 


3.84 


0.35 


2.83 


0.58 


6.71 


1.1 


3.37 


9 


cyclohexane 


5.7 


1.2 


4.5 


3.43 


0.45 


2.77 


0.44 


10.8 


1.7 


4.07 


10 


methylcyclohexane 


5.7 


1.1 


4.6 


3.74 


0.46 


2.78 


0.47 


11.0 


1.4 


3.47 


11 


benzene 


5.5 


1.1 


4.4 


3.64 


0.48 


2.77 


0.46 


12.9 


1.6 


4.90 


12 


toluene 


5.5 


1.1 


4.4 


3.77 


0.49 


2.77 


0.46 


12.8 


1.5 


4.10 


13 


ethylebenzene 


5.7 


1.2 


4.5 


3.25 


0.45 


2.76 


0.37 


10.6 


2.1 


3.60 


14 


o-xylene 


5.5 


1.2 


4.3 


3.57 


0.48 


2.77 


0.40 


12.0 


1.7 


3.73 


15 


m-xylene 


5.6 


1.2 


4.4 


3.52 


0.46 


2.77 


0.40 


11.3 


1.8 


3.54 


16 


p-xylene 


5.6 


1.1 


4.5 


3.52 


0.46 


2.77 


0.40 


11.2 


1.8 


3.51 


17 


mesitylene 


5.6 


1.1 


4.5 


3.52 


0.46 


2.77 


0.40 


11.1 


1.8 


3.13 


18 


1,2,3,4-thetra- 

methylbenzene 


5.5 


1.2 


4.3 


3.46 


0.48 


2.76 


0.40 


11.6 


1.7 


2.9 


19 


hexafluorobenzene 


5.9 


1.1 


4.8 


3.64 


0.42 


2.79 


0.42 


9.47 


1.9 


3.27 


20 


tetramethylsilane 


6.4 


1.1 


5.3 


3.77 


0.37 


2.82 


0.54 


7.15 


1.2 


2.82 


21 


diethylether 


6.2 


1.2 


5.0 


3.36 


0.39 


2.80 


0.49 


8.23 


1.5 


3.64 


22 


dioxane 


4.9 


1.0 


3.9 


4.46 


0.61 


2.76 


0.60 


19.7 


0.8 


4.07 


23 


methanol 


6.2 


1.2 


5.0 


3.21 


0.39 


2.79 


0.42 


9.37 


2.7 


8.09 


24 


ethanol 


6.0 


1.2 


4.8 


3.33 


0.41 


2.78 


0.45 


9.91 


2.1 


6.13 


25 


propanol 


5.9 


1.2 


4.7 


3.39 


0.43 


2.78 


0.45 


10.4 


1.9 


5.17 


26 


butanol 


5.7 


1.1 


4.6 


3.54 


0.45 


2.78 


0.45 


11.1 


1.7 


4.42 


27 


octanol 


5.4 


1.1 


4.3 


3.69 


0.49 


2.77 


0.45 


12.2 


1.3 


2.86 


28 


water 


4.3 


1.1 


3.2 


3.64 


0.73 


2.68 


0.49 


35.6 


1.7 


22.1 


29 


heavy water 


4.6 


1.2 


3.4 


3.33 


0.64 


2.69 


0.39 


27.9 


2.8 




30 


acetone 


5.8 


1.2 


4.6 


3.46 


0.45 


2.77 


0.47 


11.3 


1.7 


5.27 




acetonitrile 


0.0 


1 9 


41.0 




n AA 


9 77 
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il. 1 


A 


A 


32 


carbon disulfide 


5.7 


1.4 


4.3 


2.43 


0.43 


2.71 


0.35 


11.2 


2.8 


7.90 




He, 4.2 K 


18.6 


1.1 


17.5 


3.74 


0.05 


3.02 


0.98l°| 


0.094 


o.itl 


0.23 




H2, 20.3 K 


11.9 


1.2 


10.7 


3.16 


0.12 


2.95 




0.88 


5.7tf 


1.29 




N2, 77.3 K 


8.4 


1.1 


7.3 


4.05 


0.23 


2.90 


0.371°:!^ 


3.3 




3.39 




Ar, 86.4 K 


6.7 


1.0 


5.7 


4.50 


0.35 


2.86 


o.eotgil 


7.5 




4.90 



Calculating parameters of the Ps bubble we adopted that TirlcZesn is approximately equal to 0.5 
ns in all cases except He (1.9 ns), H2 (0.92 ns) and N2 (0.56 ns). 

2A/i?ws are calculated from Eq.(l.l) using respective values of a/a^ and assuming r = R {in 
accord with Eq.(4.9)). 



